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Abstract 

The partition functions of QCD2 on simple surfaces admit representations 
in terms of exponentials of the inverse coupling, that are modular transforms 
of the usual character expansions. We review the construction of such a rep- 
resentation in the case of the cylinder, and show how it leads to a formulation 
of QCD2 as a c = 1 matrix model of the Kazakov-Migdal type. The eigenval- 
ues describe the positions of N Sutherland fermions on a circle, while their 
discretized momenta label the representations in the corresponding character 
expansion. Using this language, we derive some new results: we give an al- 
ternative description of the Douglas-Kazakov phase transition on the sphere, 
and we argue that an analogous phase transition exists on the cylinder. We 
calculate the large N limit of the partition function on the cylinder with 
boundary conditions given by semicircular distributions of eigenvalues, and 
we find an explicit expression for the large N limit of the Itzykson-Zuber 
integral with the same boundary conditions. (Talk given at the "Workshop 
on high energy physics and cosmology" at Trieste, July 1993.) 
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1 Introduction 



In the last year, following the papers of Gross and Taylor |T|], there has been a 
renewed interest in the attempts to understand two dimensional QCD in the large 
N limit as a string theory. In ref. |l| it was shown that QCD2 is indeed a string 
theory by proving that the coefficients of the expansion of the QCD2 partition 
function in power series of 1/N can be interpreted in terms of mappings from a two 
dimensional surface onto a two dimensional target space. A different approach was 
pursued in refs. 0, [| f|] , where the equivalence of QCD2 on a cylinder and on a torus 
with a c = 1 matrix model was shown in various ways. In this talk we will describe 
the approach of Q, which leads to expressions for the partition functions that are 
"dual" to the usual character expansions, in the sense that they can be expanded 
in powers of the inverse coupling. This approach is inspired by the interpretation 
of the Kazakov-Migdal model as a model for high temperature QCD ||, which we 
review in Section 2. In Section 3 we will show, following ref 0, how QCD2 on 
a cylinder and on a torus is equivalent for any iV to a one-dimensional Kazakov- 
Migdal model, with the eigenvalues of the matter fields costrained to live on a circle 
rather than on a line. The generalization of these results to Yang-Mills theories 
based on other classical groups is also discussed, together with the equivalence 
of these theories with the zero coupling limit of Sutherland models |7J. In Section 
4, we derive some previously unpublished results concerning the large N limit of 
QCD2 on the sphere and on the cylinder. We show how the phase transition on 
the sphere recently discovered by Douglas and Kazakov || can be understood in 
our language, and we exhibit its close analogy with the phase transition of Gross- 
Witten and Wadia |1J. Finally, applying the techniques of ref. pT| , we argue 
that the partition function on the cylinder exhibits the same transition^. 

2 The Kazakov-Migdal model as high tempera- 
ture lattice gauge theory 

Let us consider the Kazakov-Migdal model fl3| with quadratic potential defined on 
a d- dimensional ipercubic lattice labeled by x: 

S = J2 NTr[m 2 (j) 2 (x) - J2 <P(x)U(x, x + n)4>{x + (Jt)rf(x, x + /x)], (1) 

X fJ, 

where 4>(x) is an Hermitian N x N matrix defined on the sites of the lattice, and 
U(x, x + fi) is a unitary N x N matrix, defined on the links, which plays the 
role of the gauge field as in the usual lattice discretization of Yang-Mills theories. 
The appearance of matter fields in the adjoint representation of the gauge group 
suggests the possibility that they may be thought of as the remnants of the gauge 

1 We have been informed [H that the existence of this transition has been independently 
established by D. Gross and A. Matytsin. 
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field along an extra compactified dimension, in the limit where the compactification 
radius tends to zero. Finite temperature lattice gauge theory in an example where 
this situation is realized, with time as the extra compactified dimension, and the 
temperature T as the inverse of the compactification radius. The corresponding 
Wilson action is given by 

Sw = ™E Re \ T H tr ( U n;0i) + ^ E HUn;ij) \ , (2) 
9 n [ i i<j J 

where n = (x,t) denotes the space-time position and U n] oi (Un-,ij) are the time-like 
(space-like) plaquette variables. 

In the large T limit the spacial plaquettes are of order ^ and the timelike plaque- 
ttes are frozen around U n -Qi = 1. The action (^) is invariant under a symmetry 
consisting in multiplying each timelike link by a space independent element of the 
center of the group. The order parameter of such symmetry is the Polyakov loop, 
defined by 

P(x)=trU(U Xit .,o). (3) 
t=i 

At high temperature this symmetry is broken (deconfinement phase): the time- 
like links U Xj t-o and the Polyakov loop are frozen around one of the invariant 
vacua, for instance 1, hence < P(x) >= 1 . In this situation it makes sense to 
expand P(x) around its vacuum expectation value and consider small fluctuations. 
It was shown in ref. jj that , in the large T limit, if one 

1. neglects the contributions coming from spacelike plaquettes, 

. 4>(x) 

2. expands P(x) = e i n powers of 

one finds that the finite temperature lattice gauge theory in d + 1 dimensions of eq. 
(0) coincides, up to terms of order i with the Kazakov-Migdal model of eq.(|l|), in 
d dimensions and with m? = d. 

In the large T limit the Wilson term corresponding to the spatial plaquettes 
appears as a small (order k) correction to the Kazakov-Migdal action. Some- 
thing similar was proposed in ref. |]14|, where a small Wilson term is added to 



the Kazakov-Migdal model in order to break the local Z^ symmetry and avoid 



super confinement. The results of ref. |14| are in complete agreement with the 



large temperature expansion proposed in [[|: two phases are present, according to 
whether the functional integral is dominated by configurations with small eigen- 
values of $(x) (phase 1) or large eigenvalues of (phase 2). In phase 1 the 
coefficient of the Wilson term remains small in the continuum limit and the Wilson 
action can be consistenly treated as a perturbation of the K-M action . However, 
as a theory of induced QCD this fails to explain the dominance of planar diagrams 
in the large iV limit . Instead, from the point of wiew of the large T expansion, 
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small eigenvalues of $(x) imply that < P(x) >= 1. We are in the high temperature 
deconfming phase and the spatial plaquettes are consistently depressed by the ^ 
factor. In this phase foam-like diagrams are expected to give a large contribution 
and we have no planar diagram dominance at large N. Quite the opposite happens 



in phase 2 where, according to |TJ]], the coefficient of the Wilson term is of order 



N in the continuum limit, so that the Wilson action is not any more a small per- 
turbation of the K-M term. From our point of view, if the eigenvalues of are 
large, the eigenvalues of P(x) are distributed on the unit circle, and < P(x) >= 0. 
The system is in the confined phase and T is below the critical temperature, hence 
the spatial plaquettes contribution cannot be neglected. This is the phase where 
the usual large N behaviour is expected to hold. 

In spite of this qualitative agreement with the physical picture, the naive large 
T limit outlined above fails for d > 1. We know in fact that the K-M model with 
a quadratic potential does not admit a continuum limit in that case |15||. On the 



other hand, it is known in the theory of finite temperature QCD that the original 
picture |TBJ of a complete dimensional reduction does not work due to the infrared 



divergences in the spacial plaquettes term, and one finds instead a three dimensional 
gauge theory coupled with a scalar field <ft in the adjoint representation, with a non 



trivial potential V{4>) (see for instance [17 , 18, |T 



The fundamental point for what follows is that none of these problems occur in 
d = 1: the continuum limit in the K-M model is possible with quadratic potential 
at m 2 = 1, and there are no spatial plaquettes, so that no approximation is done 
by neglecting them. Besides, we know that two dimensional QCD depends only 
on the area and the topology of the surface on which it is defined, so at finite 
temperature it should be independent of the temperature T. Hence, we expect the 
infinite temperature limit to give actually the exact result. This is shown in the 
next section by working directly on the functional integral of the continuum theory. 



3 QCD2 as a c = 1 matrix model 

We shall study first the partition function JC 2 (gi, g 2 ; t) of QCD2 on a cylinder, with 
fixed holonomies (Polyakov loops) gi and g 2 at the two boundaries (say x = and 
x = 27r). It is defined by the functional integral 

JC 2 (gi, £2; t) = J VA,VFe~ s ^6 (W(0), 9l ) 6 (W(2<ir),g 2 ) ^{gi)i>{g 2 ), (4) 

where the action S(t), as a result of the invariance of the theory under area pre- 
serving diffeomorphisms, can be written asQ 

27V r 27T 

S(t) = / dxdr Tr{F 2 - iF(d A 1 - dM - iF\A , Ax}}. (5) 

t Jo 

2 The coordinates on the cylinder are denoted by x, r. Notice that we choose the time direction 
(with coordinate r) as the compactified one. 
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Here F and A are independent fields and t = g 2 A, with A the area. W(x) denotes 
the Polyakov loop 

W{x)=Ve i S?* rMx > T \ (6) 

and S(g, h) denotes the conjugation invariant delta function on the group manifold. 
Finally ip{gi) and ip(g 2 ) are just normalization factors that depend only on the 
eigenvalues of g\ and g 2 and are chosen in such a way that the sewing of two 
cylinders corresponds to a group integration. 

In order to calculate the functional integral (§) it is convenient to choose the 
gauge OqAq = 0. In this gauge all the non zero modes of the Fourier expansion in 
t of Aq are set to zero. In 0] it is shown that the corresponding Faddeev-Popov 
determinant is exactly cancelled as a result of the integration over the non zero 
modes of A\ and F, so that in the end one obtains the following expression for )C 2 : 

K 2 {g u g 2 -t) = [vBVAe-^J?** 83 -*^ x 

x 5{W{0),g 1 )5{W{2ir),g 2 )iP(g 1 )i>(g 2 ). (7) 

where B(x) and A(x) (matrix fields on the algebra) denote the static modes of 
the A (x,t) and Ai(x,t) gauge fields respectively. Eq (|7|) is the first result of our 
analysis: the dimensional reduction of the previous section is exact and the result 
is, as expected, a KM model in one continuous dimension (the spatial dimension 
of the cylinder), where A(x) and B(x) play respectively the role of gauge field and 
matter field. There is, however, one important difference between (|7|) and the usual 
one-dimensional KM model: the boundary conditions depend on e 2mB rather than 
B. This means that the eigenvalues of B are defined on a circle rather than on a 
line or, as explained in [|], |J, that the matter field B enters the action through the 
unitary matrix defined by e 2mB . At this point, by means of standard matrix model 
techniques, the matrix B(x) can be diagonalized, the functional integral over its 
eigenvalues performed and one finally obtains 



N-l 

^^ 2;t)= ?jkwS ( " 1)Pexp 



_/V N 

27 (fa ~ 6 P{i) + 27Tl i 



i=l 



(8) 



where 9i and 4>i are the invariant angles of gi and g 2 , respectively^, while J (a) is, 
up to a phase the Vandermonde determinant for a unitary matrix, 



J (a) = Y[2sm 



Oi - Oj 



(9) 



3 For the SU (N) group the invariant angles 9i and & as well as the integers k are constrained 
by J^i ®i — J2i 4>i — J2i h = 0- These constraints are not there in the U (N) case, where however 
an extra factor (— l)^ -1 ) l ' is present at the r.h.s. 
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In (|8]) the normalization factor ip{g) has also been determined. Up to a constant, 
it is given by ip(g) = J(9)/A(9), where 9i are the invariant angles of g and A(0) is 
the usual Vandermonde determinant of a hermitian matrix. 

The formula (§) is one of the main results of our analysis. It provides for /C2 
an expansion in exponentials of 1/t, which is related to the well known character 
expansion by a non trivial modular transformation ||20|| , as explicitly checked in Q. 



By taking the 9i — > limit of eq. (ra), one obtains the modular inversion for the 



kernel on the disk, /Ci(</>; t), which has been known for many years [21 

x <pi-(i>j- Mh - h) 

{«i}=-oo Kj 



ft too iiii 2 sin 2 fa - <fc + Mh - h)} 

(10) 



x exp 



1=1 



The kernel on the cylinder, as given by eq. can be interpreted as follows: let 
9{ and 0, be respectively the initial and final positions of an ensemble of N particles 
on a circle, and let t be the time elapsed in the transition from the initial to the final 
configuration. Then Ki{.9\i 92] t) can be interpreted as the transition probability for 
such process. Notice that, because of the Vandermonde determinants, the wave 
functions are completely antisymmetric under the exchange of pairs of these parti- 
cles, which must therefore be fermions. That QCD2 would prove to be equivalent 
to a quantum theory of free fermions on a circle was already to be expected from 
eq. (0). It is known in fact that one dimensional KM model describes the singlet 



sector of an SU(N) matrix model with c = 1 [22[ [23], whose interpretation in terms 
of free fermions has been known for some time. 

The fermionic system that describes QCD2 turns out to be a particular case of a 
family of integrable one dimensional quantum mechanical models first investigated 
by Calogero p4j and Sutherland J7|. In fact, one can easily check from eq. (§) that 
£■2(91, 92i t) satisfies the differential equation 

- \r\<p) £ |^(<«) t) = o- (u) 

By taking the limit 9i — *■ in eq. ( |TTD one finds that also the kernel on a disk /Ci(</>; t) 
obeys the same equation. Upon redefinition of the kernel by JC\ — > J{(f>)K.\ = JCi, 
eq. (^) becomes the (euclidean) free Schrodinger equation for iV fermions on a 
circle. The boundary conditions are determined by this redefinition, and are just 
the boundary conditions that apply to the fermions in the zero coupling limit of 
the Sutherland integrable model related to the Lie algebra A N _i 0. 

The same equation, but with J(0) replaced by A(0), is satisfied by the corre- 
sponding quantities in one dimensional hermitian matrix models, where the place 
of K,2 is taken by the Itzykson-Zuber integral (see for instance the recent paper 
by Matytsin [11]]). Such models are then related to the zero coupling limit of the 
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Calogero integrable model pl |, in the same way as QCD2 is related to the Suther- 
land model. 

All quantities of relevance to QCD2 can readily be reinterpreted in the Suther- 
land language, as discussed in more detail in ||. For example, we already mentioned 
that K,2{0, 0; t) can be interpreted as the propagator from an initial configuration 9i 
to a final one 4>i in a time t. Similarly, the kernel on the disk /Ci(</>; t) corresponds 
to the transition amplitude from the configuration 9i — to 0j, and the partition 
function on the sphere is just the amplitude for the process in which all the fermions 
start at the origin and come back there after a time t. The meaning of the modular 
transformation relating the expression ([D to the usual character expansion is also 
clear: the integers labelling the unitary representations in the character expansion 
correspond to discretized momenta of the fermions on the circle, while eq. (|8[) gives 
the corresponding coordinate representation. 

The partition function of QCD2 on a torus can be obtained from eq. (|3|) by iden- 
tifying 8i with (pi and then integrating with the group invariant measure n« d9iJ 2 (9). 
The integral is gaussian in the invariant angles although with a complicated com- 
binatorial structure due to the sum over all permutation in eq. (|8]). This can be 
disentangled by decomposing each permutation in cycles and by calculating the in- 
tegrals corresponding to cyclic permutations. The result has a particularly simple 
form if one considers the grand canonical partition function at t = fixed. For 
SU(N) for instance one gets 

z s G T\i) = EaKV.tV 

dp if (l + qe-i^-^TY (12) 

The interpretation of this formula as the grand canonical partition function of 
a gas of free fermions with energy levels that go like n 2 for n — > oo is evident. 
The parameter (3 is the momentum corresponding to the center of mass of the 
fermions, which in SU(N) is localized (see footnote following eq. (|^)). Hence the 
corresponding momentum {3 is completely undetermined. 

All the results obtained in this section can be generalized to the case of an arbi- 
trary simple gauge group ||. One obtains that the matrix model describing a two 
dimensional Yang-Mills theory on a cylinder is a KM model where the matter fields 
are in the fundamental representation of the gauge group. Since the string interpre- 
tation of this matrix model requires expanding the group matrices as exponentials 
of the matrices on the algebra, one has an orientable string theory for the SU(N) 
group, while if the gauge group is SO(N) or Sp(2N) the worldsheet of the string 
may be both orientable or nonorientable 
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4 Large N phase transitions on the sphere and 
on the cylinder 



In this section we report on some new results concerning mainly the large N limit 
of the kernel on the cylinder and the partition function on the sphere, where a large 
iV phase transition was recently found by Douglas and Kazakov ||. Unlike the 
authors of f|, we work in configuration space and consider the expression for the 
partition function on a sphere of area A, obtained by sewing disks of areas t and 
A-t, 



z g=q (n,a) = I* d9 i j 2 {e))c 1 {e,t)K: 1 {e,A-t) 

Jo 

N 2 - l 

't(A-t)\ 2 




d\i J2 A(A)A(A + 2t™) x 

{rii} 



Af (A, + 2 7 rn i ) 2 
2t 2{A - 1) 



(13) 



where we used eq. ( pTOD for /C^ and, in the first kernel, we replaced 6^ + 2irmi 
with a variable Aj defined on the whole real axis. We want to find the eigenvalue 
distribution that, for any given t, dominates the integral in eq. ([13]) in the large N 
limit. The crucial point is that, if the eigenvalue distribution is, for any t, confined 
in the interval (— tt,tt), then we are entitled to consider in eq. flTB]) only the term 
corresponding to n« = 0, for all i. The other terms in fact describe the winding of 
the eigenvalues, which is possible, at iV = oo, only if the eigenvalue distribution 
spreads at some t over the whole circle. The term n« = describes a gaussian 
matrix model whose eigenvalue distribution is given by Wigner's semicircle law 



p(A,r) = — -Vr 2 - A 2 



where the radius r of the distribution is given in terms of t and A by 



lt(A-t) 



A 



(14) 



(15) 



The maximum value of the radius of the eigenvalue distribution occurs at t = =j, 



and it is given by r max = y A. According to the previous argument, r max must not 
exceed tt, which implies A < ti 2 . In complete agreement with the results of ref. || 
we find then two phases: one for A < tt 2 , where the partition function on the sphere 
is completely described in the large TV limit by a gaussian matrix model, and the 
winding modes can be neglected. In this region the Sutherland and Calogero models 
are indistinguishable at large N. The other phase, for A > tt 2 , is characterized in 
configuration space by winding numbers different from zero. It is remarkable that, 
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in configuration space, the mechanism of the phase transition discovered by Douglas 
and Kazakov is the same as the one of the Gross- Witten-Wadia phase transition 



1C 



This reasoning can be extended to the cylinder if one realizes that, from the 
point of view of the dynamics of the eigenvalues, the sphere is just a cylinder with 
special boundary conditions, namely that all eigenvalues must be concentrated at 
the origin at the initial and final times. To study the cylinder in the large N limit, 



in the gaussian phase, one needs to generalize eq. (15), to find the time evolution 
of the radius of the semicircular eigenvalue distribution for arbitrary boundary 
conditions (that is allowing for non-zero initial and final radii). This can be done 
by making use of the results of a recent paper by Matytsin |TTJ , which studies the 
time evolution of the density of eigenvalues in a system satisfying eq. (|TTD , with J(0) 
replaced by A(0). The result is that the evolution is determined by the Das-Jevicki 
hamiltonian pEfl, an d the corresponding equations of motion can be written as the 



Hopf equation 

f + ^°- < ie > 

where f(X,t) = v(X,t) + z7rp(A,t), and p(X,t) is the density of eigenvalues Aj at 
time t. According to our discussion, the replacement of J(4>) with A(0) in eq. ([TT| ) 
does not affect the large N limit in the gaussian phase. Therefore, the semicircular 
distribution (|HD, with radius given by (|i~5|), must be a solution of ([R]). Indeed, if 
one inserts in (|16D the ansatz (|T4]) , with r = r(t), one finds that flUf ) is solved if 



r + - = 0. (17) 



The general solution of this differential equation is easily found, and it generalizes 
naturally eq. (|15"D, as 



r(t) 



! (t + a)((3-t) 



V a + 13 

Given arbitrary boundary conditions r(0) = 70 and r(A) = 71, one can find from 
(|j~8|) the corresponding values of a and j3. In particular, the sphere corresponds to 
7o = 7i = 0, which implies a = and (3 = A, in agreement with ([151) . Eq. (|14]) , 
with r given by (|T8|), gives then the classical trajectory for the density of eigenvalues 
that dominates, in the large iV limit, the functional integral corresponding to the 
kernel /C2, with area A and boundary conditions given by Wigner's distributions of 
radii 70 and 71. The actual value of /C2 in the large N limit is given by 

/C 2 ( 7o ,7i;^) = e- iv25 — , (19) 
where Sdass is the action calculated on the classical trajectory, 

S class = lJ Q A dt(f* + ±) (20) 
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with r given by (|I%D. 

Now it is easy to see that the same mechanism that leads to a phase transition 
on the sphere operates also on the cylinder. The solution ([HJ) is the large N limit of 
the QCD2 kernel on a cylinder only if r(t) < tt for any t on the trajectory. For any 
given boundary conditions (that is for any given 70 and 71), the maximum value 
of r(i) increases as the area A increases (this is because f < from (|I7]) ), so it is 
bound to hit the value tt at some critical value of A, which can easily be determined 
from the initial conditions. There a phase transition occurs, in analogy with case 
of the sphere. 

Our solution for the large N limit of /C2 of course applies also to matrix models 
where the eigenvalues live on the real axis rather than on a circle, but in that case 
with no restrictions on the values of r(t), and no phase transitions. In other words, 
our solution applies without restrictions if one is dealing with a Calogero rather than 
a Sutherland model. In this case the solution flT9|) is equivalent to the calculation 
of the Itzykson-Zuber integral in the large N limit |pT[ , with Wigner's distribution 
of the eigenvalues. We find 

-^lnJ 7Z (7o,7i) = 

/ / 2 2" \ 

+ 0(1/N), (21) 

where I/z(7o,7i) is the Itzykson-Zuber integral, calculated with respect to two sets 
of eigenvalues that in the large N limit have Wigner distributions of radii 70 and 




71. Notice that, for 70 = 71, the expression (|2l|) reduces to the one calculated by 
Gross in [ |15| . 

We see that the "dual" expansions for the partition functions of QCD2, which 
have been determined so far for the disk, the sphere, the cylinder and the torus, are 
valuable tools for the study of the theory, both for finite N and in the large N limit. 
They provide a physical insight which is complementary to the one given by the 
character expansions, just as configuration space is complementary to momentum 
space, and they have a very direct interpretation in terms of Sutherland fermions, 
which leads to a simple understanding of the phase transitions on the sphere and 
on the cylinder. 
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